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ABSTRACT 

We report on time-resolved optical imaging of the X-ray binary SAX J1808.4— 3658 during its 
quiescent state and 2008 outburst. The binary, containing an accretion-powered millisecond pulsar, 
has a large sinusoidal-like modulation in its quiescent optical emission. We employ a Markov chain 
Monte Carlo technique to fit our multi-band light curve data in quiescence with an irradiated star 
model, and derive a tight constraint of 43°±3° on the inclination angle of the binary system. The pulsar 
and its companion are constrained to have masses of 1.6±°; 3 * M Q and 0.07±°;°i M , respectively. We 
find that the accretion disk had nearly constant fluxes over a ^500 day period in the quiescent state 
our data covered, but was caught to have started brightening 1.5 months before the 2008 outburst. A 
change of the disk to a steeper spectrum was also seen at the onset of the brightening. Variations in 
modulation during the outburst were detected in our four observations made during 7-12 days after 
the start of the outburst, and a sinusoidal-like modulation with 0.2 mag amplitude changed to be flat 
with a small amplitude of 0.1 mag. The modulation variations are discussed. We estimate the albedo 
of the companion during its quiescence and the outburst, which was approximately and 0.9 in the 
former and latter states, respectively. This large difference probably provides additional evidence for 
that the neutron star in the binary turns on as a radio pulsar in quiescence. 

Subject headings: binaries: close — stars: individual (SAX J1808.4— 3658) — X-rays: binary — stars: 
low-mass — stars: neutron 



1. INTRODUCTION 

Among ^200 known low-mass X-ray binaries 
(LMXBs), SAX J1808.4-3658 (hereafter J1808.4) 
stands out. Not only was it the first discovered 
accretion-powered millisecond pulsar (APMP) system 
( Wijnands fe van der Klisl 119981 ). but also this transient 
system is relatively bright during quiescence and has 
relatively frequent outbursts, allowing detailed s tudies 
of its various properties (see lHartman et al.l 120081 and 
references therein). One possible property of the system 
that is extremely interesting and thus far unique amo ng 
over a dozen known APMP systems (jPatrunol 120101 ) is 
that the millisecond pulsar switches its energy source 
from accretion during outbursts to rotation during qui- 
escence. This possible property is indicated by bright, 
large-amplitude optical modulation seen in the binary 
during quiescence, which can not be caused by X-ray 
heating of the inner face of the companion star because 
the quiescent X-ray luminosity L\ (~ 5 x 10 31 erg s _1 ; 
ICampana et al.1 [20021 : iHeinke et al.l I2007T ) is two orders 
of m agnitude lower than that required (|Burderi et al.l 
2003). Considering that the neutron star in the binary 
turns on as a radio pulsar in quiescence, its rotational 
energy output (so-c alled spin-down lumi nosity) L s d is 
~ 9 x 10 33 erg s _1 ([Hartman et al.ll2008D . This energy 
output would presumably be in the form of a pulsar 
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wind and sufficient to illuminate the inner face of the 
companion star. 

To fu lly examine the po ssibility, Deloye et al. (2008; 
see also IHeinke et a l. 2009) have observed J1808.4 simul- 
taneously at X-ray and optical energies, and from the ob- 
servations they have confirmed the inconsistency between 
the large-amplitude optical modulation and low X-ray lu- 
minosity. Applyin g a phase-coherent timing technique, 
I Wang et al.l ()2009l ) have accurately determined the pe- 
riod and phase of the optical modulation and concluded 
that the optical periodicity is highly consistent with the 
X-ray orbital epheme ris (derived from pulsar timing; see 
lHartman et al.l [20081) . The results have excluded other 
possible origins (such as the accretion disk) and estab- 
lished that the modulation does arise from the compan- 
ion star. 

Following these studies, it should be interesting to com- 
pare the modulations in an outburst and quiescence since 
presumably the former is caused by X-ray heating, and 
the latter by pulsar wind heating. W hen J1808.4 was re- 
ported to have its expected outburst (|Gal lowav 2008) at 
the end of 2008 September (jMarkwardt fe Swankl T2008) . 
we organized optical observations of the source and ob- 
tained its light curves during the outburst. In this paper, 
we report on t he results from the o bservations. 

In addition, iDelove et al.1 ([20081 ) have shown that by 
fitting the quiescent optical modulation in J1808.4 with 
an advanced binary light curve model, in which lights 
from an accretion disk and an irradiated donor are 
considered, important constraints on the binary system 
(such as the inclination and neutron star's mass) can be 
derived. It has also been realized by them that in order 
to have better constraints, simultaneous multi-band light 
curves are needed. Using data recently obtained as well 
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as tho se reported in lDelove et al.l (|2008l ) and I Wang et al.l 
we have conducted fitting to the modulation thor- 
oughly. We used a binary light cur ve code, adapted from 
that developed by van Kerkwijk (iStappers et alJl200lD. 
that is similar to the ELC program (jOrosz fc Hauschildtl 
12000( 1 used bv iDelove eFall (|2008| ). We also employed 
a Markov chain Monte Carlo (MCMC) technique that 
allows for simultaneously fitting the multi-dimensional 
light curve model to the data and identifies the best-fit 
parameters as well as their confidence intervals. 

We describe our time-resolved photometry of J1808.4 
and data reduction in § [2j Details of our binary light 
curve model and MCMC fitting of the obtained light 
curves, and fit results are given in § [3j We discuss inter- 
esting results from our observations in § [4j and summa- 
rize in § [5] 

2. OBSERVATIONS AND DATA REDUCTION 

2.1. Gemini Imaging 

Time-resolved imaging of J1808.4 was carried out with 
Gemini South Telescope on 2008 August 4 and 7. To ob- 
tain nearly simultaneous multi-band light curves, Sloan 
g' and i' filters were alternately used for imaging in each 
night. The detecto r was the Gemini Multi-Object Spec- 
trograph (GMOS; iHook et al.l 120041 ) . which consists of 
three 2048x4608 EEV CCDs. In our observations, only 
a section of 300x300 pixel 2 in the middle CCD (CCD 
02) was used. The pixel scale is 0.073"/pixel. 
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Fig. 1. — Gemini South g' (upper panels) and i' (bottom panels) 
light curves of J1808.4. Light curves of a comparison star (trian- 
gles) are also plotted. 

The total observation time in the first night was 3 
hours, with 23 and 24 useful g' and i' images respec- 
tively made. The exposure times were 200.5 sec in g' 
and 150.5 sec in i' . The observing conditions were vari- 
able, with the seeing [full- width half-maximum (FWHM) 



of the point spread function (PSF) of the images] increas- 
ing from 0.6" to 0.9" over the course of the observation. 
On August 7, the total observation time was 5.5 hours, 
during which 47 and 48 useful g' and i' images, respec- 
tively, were made. The exposure time at each band was 
the same as that in the first night. The observing condi- 
tions were relatively stable, with the average seeing being 
-0.7". 

2.2. Canada France Hawaii Telescope ( CFHT) Imaging 

The sta rting of the 2008 outburst o f J1808.4 was re- 
ported bv lMarkwardt fc Swankl (|2008D on 2008 Septem- 
ber 22. We subsequently requested CFHT Target-of- 
Opportunity observations of the source. The observa- 
tions were carried out in the queued service observ- 
ing mode on 2008 September 29, 30, October 3, and 
4. The detector was the wide-field imager MegaCam, 
which consists of 36 2048x4612 pixel 2 CCDs. The field 
of our target was imaged on CCD 22. The pixel scale is 
0.187"/pixel. A Sloan r' filter was used for imaging. 

During the four nights, we obtained 60, 63, 84, and 72 
images, respectively. The exposure time for each image 
was 20 sec. The total observation times were 1.1, 1.2, 
1.6, and 1.4 hours, respectively. The observing condi- 
tions in the first night were good, with an average seeing 
of 0.8". In the second night, the seeing condition was not 
good, having a large range of 0.8-1.8" and an average of 
1.1". Therefore 6 images with large seeing values were 
excluded from the data. The third night had good ob- 
serving conditions, with an average seeing of 0.7". The 
observing conditions in the fourth night were relatively 
variable, with the seeing in a range of 0.6-1.2"and an 
average of 0.8". 

2.3. Data Reduction and Photometry 

We used the IRAF packages for data reduction. The 
images were bias subtracted and flat fielded. We per- 
formed PSF-fitting photometry to measure the bright- 
nesses of the source a nd other in-field stars . A photom- 
etry program D0PH0T (jSchechter et al.lll993h was used. 

Before photometry, we positionally calibrated our 
Gemini images (made during the quiescent state) to re- 
duce possible contamination from the two nearby stars, 
since the two stars had similar bright nesses and were 0.6 " 
and 1.0" away from the target (e.g.. iWang et al.l f2009). 
We first made a reference image by combining three best- 
quality i' images on August 7. All Gemini images were 
then calibrated to this reference image. We determined 
the positions of our target and the two nearby stars in the 
reference image and fixed them at these positions for pho- 
tometry of the Gemini images. Differential photometry 
was performed to eliminate systematic flux variations in 
the images. Six isolated, non- variable bright stars in the 
field were used. The brightnesses of our targets and other 
stars in each frame were calculated relative to the total 
cou nts of the six stars. Standard stars 95-100 and 95- 
96 (ILandoltlll992T ) were observed on August 4 and used 
for q' and i' fl ux calibration, respectively. To convert 
ILandoltl (11991 Ve ga magnitudes of the standard stars to 
Sl oan filter magni t udes, transformation equations given 
bv iFukugita et al.l (|1996| ) were used. 

For CFHT data, because J1808.4 was in outburst, ap- 
proximately 30 times brighter than in quiescence, and 
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Fig. 2. — CFHT r' light curves of J1808.4 in 2008 outburst (diamonds), with the last two up-shifted by 0.3 mag. For comparison, light 
curves of a comparison star are plotted as triangles. The X-ray ephemeris of the pulsar is shown as dotted curves, and positions of superior 
conjunction of the companion are indicated in the last two panels (dash dot lines). Our best-fit model light curves are shown as dashed 
curves. 



the exposure time was short, no effort was made to sep- 
arate the two nearby stars from our target. Differential 
photometry was also performed, with 10 in-field bright 
stars used for calibrating out systematic flux variations 
among the images. No standard stars were requested in 
our CFHT program. Using 8 in-field bright stars , we flu x 
calibrated our r' images to those in Wangj5t al. (2009). 

2.4. Light Curve Results 

The obtained Gemini and CFHT light curves are 
shown in Figure [T] and Figure [2] respectively. We note 
that during our Gemini observations, the source was 
brightening. Over two days, the source's g' and i' bright- 
nesses changed by 1.3 and 0.9 mag, respectively. In addi- 
tion, the amplitudes of modulation decreased from ~0.4 
mag to ~0.2 mag. 

At the time of the outburst, J1808.4 was visible to 
CFHT for only approximately one hour. Therefore each 
CFHT light curve covered half of the binary orbit. It can 
be seen that the first two light curves are sinusoidal-like 
with an amplitude of 0.2 mag, similar to those seen in 
quiescence. The latter two light curves however changed 
significantly to be flat (or more-or-less triangle-like) with 
an amplitude of 0.1 mag, and the average source bright- 
ness at the time decreased by approximately 0.3 mag. 

3. FITTING 

3.1. Light curve model 

We fitted the data using a light curve model for 
irradiated companions in binaries. The model was 
adapted from the co de developed by van Kerkwijk (e.g. 



Stappcrs et al. 



20011) . and has been recent ly used for two 



other systems ([van Kerkwijk et ~aT] |2010bllah . The model 
is sim ilar to the ELC pr ogram (lOrosz fe HauschildtJ 
I2000D that was employed bv lDelove et all (|2008l) for their 
earlier work on J1808.4. 



In this model, the surface of the companion is con- 
structed by solving the hydrostatic equilibrium equation 
for a rotating body, which in our case was tidally locked 
to the neutron star. We implemented a new stellar grid 
parametrization that uses a triangle tessellation obtained 
from the subdivision of the primitives of an isocahedron. 
This parametrization yields uniform coverage on the pro- 
jected sphere surrounding the star and allows for better 
performances of the code. We scaled the star so that 
its size corresponds to that of the Roche lobe, which is 
justified from the fact that J1808.4 regularly undergoes 
outbursts involving mass transfer from the companion. A 
base temperature for the companion is chosen and effects 
of gr avity darkening accounted using Lucy's law ()Lucvl 
1967), with a power-law index of 0.08. To this, we ap- 
ply the irradi ative flux from the neutron star in the way 
prescribed bv lOrosz" fe Hauschi ldt (2000). Finally, we in- 
tegrated over the visible surface the emerging flux in a 
selected photometric filter, w hich was interpolated f rom 
NextGen model atmospheres (lHauschildt et al]|1999f) , in 
order to obtain the total flux at a given orbital phase and 
for specific orbital parameters. 

We approximated the disk contribution by adding an 
orbital-independent, constant value to the companion's 
flux. This approach is well-suited given that what we 
obtained were the long-term changes of the disk contri- 
bution. The specifics of the spectral energy distribution 
relating the different photometric bands can be inferred 
after the fact. We also made the implicit presumption 
that the system is seen sufficiently face-on so that there 
is no mutual eclipse/shadow arising between the disk and 
the companion. This presumption is motivated by the 
fact that light curves of J1808.4 in quiescence are sym- 
metrical and do not show obvious signs of flux changes 
that could result from the interaction between the two 
components. Moreover, the estimated orbital inclination 
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from previous work (e.g. lDelove et all I2008T ) as well as 
that inferred from this work (see £ )3.5I) indicates that the 
companion is unlikely to be obscured by the disk for any 
reasonable disk scale height. 

3.2. Bayesian inference 

We used a Bayesian framework to make statistical in- 
ferences on our light curve modeling parameters. The 
Bayesian inference on a set of parameters 9 is obtained 
from their posterior probability, which is conditional to 
the prior knowledge, i, that we have of these parameters 
and the experimental data D. The posterior probability 
can be written as the product of the priors times the like- 
lihood of the data: p(6\I,D) = p{6\I) p(D\9, I)/p(D\I), 
where the denominator is a normalization factor to en- 
sure that the posterior prob ability integrate d over the 
parameter space is unity [see iGregorvl ((2005) for an in- 
troduction to Bayesian analysis]. 

3.3. Model parameters 

INCLINATION AND RADIAL VELOCITY 
SEMI-AMPLITUDE 

The X-ray timing of J1808.4 already provides accu- 
rate measurements of its orbital period and the projected 
semi-major axis of the neutron star. In order to fully 
parametrize the system (e.g. masses, separation, etc.), 
two additional quantities are required. For convenience 
reasons, we chose the orbital inclination, i, and the pro- 
jected radial velocity semi-amplitude of the companion 
-P^Comp, because of their close connection to observables. 
For the orbital inclination, we chose the priors to be flat 
in cosi. 

The use of spectroscopic observations to constrain the 
companion's radial velocity is made quite difficult due 
to its faintness in quiescence. On the other hand, while 
the system gets much brighter during outbursts, its emis- 
sion mainly contains features arisi ng from the ac c retion 
disk. ICornelisse et all (|2009fl and lElebert et all (|2009f) 
respectively tracked the companion's motion using spec- 
tra of J1808.4 obtained in the 2008 outburst. Here we 
used results given by the latter, because of large uncer- 
tainties in the results in the former, due to limited data 
and poor S/N of the data. Through Doppler tomog- 
raphy of the N III A4640.64 emission line associated to 
the Bowen blend, a fluorescence feature caused by the 
heating of the com panion by energetic UV photons orig- 
inating in the disk, Elcb ert et al.l (|2009f) found ifeowcn = 
324 ± 15kms _1 . However, ifsowen follows the average 
location of the Bowen emission, which is essentially pro- 
duced on the irradiated hemisphere of the star and has a 
different veloc i ty fro m that of the center of mass, -Kcomp- 
lElebert et al.l ()2009D discussed the "K-correction" factor 
required to convert the Bowen velocity to the center of 
mass and estim ated that -ffpoma = 370±40 kms~ 1 , based 
on modeling by iMunoz-Da rias et al.l (|2005l ). 

We included the if Bowen constraint as a Gaussian prior 
in our Bayesian inference. For every realization of the 
model, we determined the velocity of the companion's 
light center by calculating the average velocity over the 
visible surface, weighted by the flux contribution of each 
surface element. Since the irradiated side of the com- 
panion is much hotter than that of the back side (6000 K 
vs. 3000 K, respectively), the flux contribution from 



TABLE 1 

Fitted and Derived Model Parameters 



Parameter 


Median 


68.3% interval 


99.7% interval 


Inclination (degree) 


43 


40,46 


34,53 


Kn f'kms -1 ) 
J1 Oomp V o ) 


378 


372 384 


359 397 


Distance modulus 


12.78 


12.73, 12.84 


12.63, 12.94 


Absorption (J-band) 


0.21 


0.19, 0.23 


0.15, 0.27 


MJD 54167.35 








Day-side Temp. (K) 


7290 


7120, 7460 


6800, 7830 


Disk g' (pjy) 


3.2 


2.6, 3.8 


0.9, 6.2 


Disk i' (fi3y) 


6.9 


5.8, 8.1 


2.2, 12.4 


MJD 54169.35 








Day-side Temp. (K) 


7010 


6850, 7170 


6550, 7500 


Disk g (/ijyj 


3.7 


3.1, 4.4 


1.2, 7.1 


T)it;lf 9 f ( II 1\A 
-L/ioJv i \Ltoy ) 


8 


6.7 9.3 


2.7 14.1 


MJD 54597.32 








Day-side lemp. (K) 


6780 


6600, 6960 


/ ' < *W ' i 1 Tl /lit 

6260, 7340 


r>i<ik T 1 ( ii Tv^t 


5 1 


4.3 6.0 


1.7 9.3 


MJD 54598.34 








Day-side Temp. (K) 


6690 


6520, 6870 


6190, 7230 


Disk r' (fj.3y) 


5.3 


4.5, 6.2 


1.8, 9.6 


MJD 54601.38 








Day-side Temp. (K) 


6700 


6530, 6880 


6210, 7250 


Disk r' (/^Jy) 


5.8 


4.9, 6.8 


2.0, 10.5 


MJD 54682.08 








Day-side Temp. (K) 


6600 


6440, 6770 


6140, 7120 


Disk g' (fijy) 


4.6 


3.8, 5.4 


1.6, 8.7 


Disk i' (/ijy) 


6.6 


5.6, 7.7 


2.3, 11.6 


MJD 54685.12 








Day-side Temp. (K) 


6520 


6320, 6720 


5960, 7120 


Disk g' (pjy) 


16.6 


14.0, 19.6 


6.4, 31.2 


Disk i' (>Jy) 


15.5 


13.2, 18.0 


5.9, 26.9 


Outburst data 


MJD 54738.21 








Day-side Temp. (K) 


14860 


14250, 15460 


12980, 16620 


Disk r' (^tjy) 


130.7 


129.9, 131.6 


126.0, 135.6 


MJD 54739.21 








Day-side Temp. (K) 


17540 


17050, 18040 


16150, 18950 


Disk r' (/ijy) 


129.2 


128.3, 130.1 


124.0, 133.6 


MJD 54742.21 








Day-side Temp. (K) 


9110 


8780, 9440 


8180, 10060 


Disk r' (^tJy) 


108.0 


107.3, 108.7 


103.7, 111.4 


MJD 54743.21 








Disk r' (AtJy) 


106.9 


106.0, 107.7 


101.8, 111.0 


Inferred Parameters 


Mass Comp. (M ) 


0.07 


0.06, 0.08 


0.04, 0.12 


Mass NS (M ) 


1.60 


1.35, 1.91 


0.99, 2.88 



the non-irradiated side in the (/-band becomes negligible 
and the inferred velocity should be a good approxima- 
tion o f if Bowen- As expected from iMuhoz-Darias et al.l 
(2005), we found that ifcomp ~ 370 km s -1 is required 
to satisfy the velocity constraint. Such value is also 
consistent with the semi-analytical K- c orrecti on rela- 
tionship derived by Ivan Kerkwiik et~aTI ()2010al ) for the 
companion to PSR B1957+20, which shares similarities 
to J1808.4. From analytical and numerical modeling 
of this irradiated companion, they determined that the 
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ratio of the observed projected light-center velocity to 
the projected center-of-mass velocity follows K Q \ 1S /K2 = 
1 — fcsRnosc/o-2, where Rnosc/0,2 is the radius of the nose 
of the companion in units of separation to the system's 
barycenter, and / c g ~ 0.6, a parameter relatively inde- 
pendent of the system and companion properties such 
of orbital inclination, filling factor, and surface tem- 
perature as long as the companion is strongly irradi- 
ated and nearly fills its Roche lobe. Hence, for J1808.4 
■Rnose/^2 ~ 0.2 which, again, implies a measurable veloc- 
ity comparable to the above value. 

BASE AND DAY-SIDE TEMPERATURES 

The base temperature of the companion was fixed 
using the tem perature-mass rela tionship for low-mass 
brown dwarfs (|Delove et al.l 120081 ). The temperature, as 
a function of mass alone, ranges between 2000-3000 K. 
From the photometric colors near the inferior conjunc- 
tion of the companion, we confirmed that the tempera- 
ture is consistent with the range. In addition, keeping 
it a free parameter did not improve the light curve fits 
significantly. 

We used flat priors for the companion's day-side tem- 
perature since its surface albedo is uncertain as well as 
the amount of incident energy approximated by the pul- 
sar's spin-down energy. Although it is unlikely, we tested 
modeling the whole dataset with a constant day-side tem- 
perature. We found such setting did not result in good 
fits. A separate day-side temperature for each observa- 
tion was set in order to account for variability. 

DISTANCE MODULUS AND REDDENING 



iGallowav &: Cummingl (|2006| ) obtained a distance to 
JI808.4 using the modeling of the X-ray flux during the 
2002 October outburst and found d = 3.4 — 3.6 kpc. We 
took this measurement to set a Gaussian prior on the 
distance modulus DM = 12.72 ± 0.15. Here we chose a 
larger, more conservative uncertainty in order to compen- 
sate for possible systematic uncertainties in the distance 
estimate. 

For the reddening, we also chose a Gaussian prior and 
set the J band abs orption to A j = 0. 21 ± 0.05 using the 
value derived from lWang et al.l ()2009f ). 

DISK FLUX 

The same as for the day-side temperature, we used an 
independent disk flux value with a flat prior for each set 
of continuous observation, as well as each photometric 
band. We designed our light curve fitting algorithm to 
optimize the disk flux value every time the companion 
parameters are changed. Since the disk contribution is 
an additive, linear problem in the flux system — and 
therefore a very simple non-linear problem in the mag- 
nitude system — solving for the best disk flux comes to 
the cost of only a few additional CPU cycles compared to 
the lengthier optimization of the companion parameters. 

In the 2008 August 7 observation (Figure [lj , a notice- 
able rise in the overall flux level is visible and likely indi- 
cates that J1808.4 was coming out of quiescence. Clearly, 
a constant flux value cannot account for this change and 
we allowed for a small linear variation of the disk contri- 
bution over the course of the observation. 



3.4. MCMC fitting 

The light curve fitting is a multi-dimensional nonlinear 
optimiza tion problem tha t we tackled using an MCMC 
method (Gilk s~et al.lll995l ). This is particularly efficient 
to make Bayesian inference since the full posterior dis- 
tribution is evaluated. In the limit that the observed 
magnitude errors are normally distributed, the likelihood 
function can be written as p(D\0I) = exp(— x 2 /2), where 
X 2 is the conventional chi-square. 
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Fig. 3. — Folded quiescent light curves, with g' , r' , and %' data 
shown in the top, middle, and bottom panel, respectively. Fol- 
lowing the data sequence given in Table [T] the data are plotted as 
diamonds, triangles, squares, and plus signs, respectively, while the 
best-fit model light curves are plotted as dotted, dashed, dash-dot, 
and long dash curves, respectively. Two cycles are displayed for 
clarity. 

For the MCMC, we chose the stretch move algo- 
rithm that is part of a family of MCMC methods called 
ensemble samplers with a ffine invariance described by 
iGoodman fc Weard fj2009f ) . In a nutshell, the algorithm 
is executed by simultaneously running several chains, all 
initialized at random locations of the parameter phase 
space. Every step of the MCMC, a move is proposed for 
each chain by choosing a complementary chain at random 
and drawing point along a line passing through the last 
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Fig. 4. — One- and two-dimensional distributions of the key parameter values resulting from our MCMC fitting. Values of the median, 
mean, and standard deviation of each parameter are given. The masses are inferred values from the neutron star mass function, the orbital 
inclination, and the projected velocity. 



position recorded in the two chains. How big a step away 
from the previous position is determined using a distri- 
bution which is affine invariant [i.e., g(z~ 1 ) = zg(z)]. 
The proposed move is accepted with a probability that 
is slightly different than that of the usual Metropolis al- 
gorithm i n order to satisfy th e deta iled balance require- 
ment [see lGoodman fc W care (2009) for the details]. An 
advantage of the stretch move algorithm is that the ac- 
ceptance rate, which controls the efficiency of the chain, 
can be tuned using a single parameter in the proposal 
distribution, as opposed to one per dimension for a ran- 
dom walk Metropolis algorithm. It should also perform 
better for highly skewed and badly scaled distributions. 

We ran the stretch move algorithm using 18 simultane- 
ous chains for a total of 100 000 steps each. The proposal 
distribution's tuning parameter was determined from a 
previous trial run and the obtained acceptance rates of 
the chains were found to be in the range 20 — 30%, as pre- 
scribed bv lRoberts fc Ros"en thal (200j]). For subsequent 
analysis, we discarded the first 10 000 points (e.g., the 
burn-in period) and thinned the remaining by keeping 
every other 30 points in order to reduce the autocorrela- 
tion. In order to ensure that our MCMC has converged 
we visually inspected the t race of the p arameters and 
performed the Ceweke test (|Gewekdll992D by comparing 
the mean of 10%-subsections of the first half of the chain 
to the remaining second half. 

3.5. Results 

We originally aimed to compare the modulations in 
quiescence and outburst by fitting the light curves sepa- 
rately. However the combination of the short light curves 
during the outburst and their relatively large uncertain- 



ties does not provide tight constraints. While we also 
tested to fit all data together, for the same reason the 
fit results were hardly changed compared to those from 
fitting the quiescent data along. We therefore first fit the 
quiescent data, and the reduced \ 2 f° r our best-fits are in 
a range of 0.3-4. In our case, x 2 values reflect goodness- 
of-fit relatively since smaller values are obtained from 
fitting data with poorer quality (i.e., larger uncertain- 
ties). In addition, it has been known from previous stud- 
ies that significan t intrinsic variations appear in modu - 
lation of J1808.4 (jDelove et al.ll2008t IWarig et alJl2009f ). 
In Figure El we show all folded quiescent light curves as 
well as the best-fit models. Table Q] presents a summary 
of the fit results as well as several additional quantities 
inferred from the marginalized posterior probabilities of 
our irradiated companion modeling. The one and two- 
dimensional posteriors are displayed in Figure |4] 

Our inferred orbital inclination, i = 43° ± 3°, i s con - 
siste nt with previ o us est imates from iWang et al.l (|200lD 
and iDelove et al.l (|2008D (though on the lower side of 
the former) and improved substantially despite the day- 
to-day disk and irradiated-side temperature variability, 
which probably impair our ability to combine multi- 
epoch data obtained in different bands. Consequently 
we constrained the neutron star to have a mass of 
1.60±^q'25 ^0 an d the companion star to have a mass 
of 0.07±0.01 M©. The strong constraint on the neutron 
star's mass comes from the K comp measurement, whic h 
prefers a non- heavy neutron star ( [Elebert et al.l F2009) . 
Our fitting has narrowed the possible mass ranges for 
both the neutron star and companion. 

The results of the distance modulus and extinction 
are naturally in agreement with previous estimates since 
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they were used as priors for our inference. Neverthe- 
less, we can assess that the new data has an impact at 
least comparable to that of the priors since the posteriors, 
which are very close to being normally distributed, have 
standard deviations only 33% and 40%, respectively, as 
large as that of the priors. 

Similarly, our inferred companion's projected radial 
velocity semi-amplitude, 378 ± 6kms . is in agree- 
ment with the value, 370 ±40 km s -1 , from lElebert et ahl 
(2009). We observed that uncertainty on i^comp, is only 
15% as large as that of the priors, which implies that 
the data significantly contribute to constraining the com- 
panion's velocity. The most likely reason is that given 
the distance and absorption relatively well-determined, 
and the fact that the companion fills its Roche lobe, the 
observed luminosity constrains the actual radius of the 
companion which is tied up to its orbital velocity. 

We found that the day-side temperature of the com- 
panion appears to decrease over time before the 2008 
outburst (see Figure [5]). A linear model with a slope 
of -1.06 K day" 1 and an intersect of 7149 K at MJD 
54167.3503 is strongly favored over a constant temper- 
ature model. This behavior is relatively surprising as 
one might have naively expected more activity since the 
system was going towards an outburst. Considering that 
the companion was still irradiated by the pulsar wind, a 
possible explanation is that the accretion disk was thick- 
ening over the time, thus gradually blocking more part of 
the pulsar wind. On the other hand, the behavior of the 
disk flux contribution (Figure [6]) is consistent with be- 
ing constant in each band over time, except for the last 
observation in which the flux increased by a factor > 2 
within 3 days. Detailed analysis of the disk components 
is given below in 
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Fig. 5. — Day-side temperature of the companion star in J1808.4 
as a function of time. The temperature shows a trend of decreas- 
ing, which can be described by a linear model having a slope of 
-1.06 K day" 1 . 

3.5.1. Outburst results 

Because of the reasons given above, to study the out- 
burst light curves, we fit them by fixing all parameters of 
the binary system at the best-fit values obtained from our 
fitting to the quiescent data. The free parameters were 
only the day-side temperature of the companion and disk 
flux. In addition, to have sufficient constraints, the last 
two sets of the data were assumed to have one day-side 



temperature value (the average values of the data do not 
have a significant change). Our model could generally fit 
the light curves, as the reduced x 2 are hi a range of 1-2, 
similar to those obtained from fitting the quiescent data. 
From fitting, we found that the day-side temperature was 
in a range of 9000-17000 K and disk flux approximately 
20 times larger than those in quiescence. These fit results 
are given in Table [1] A test to use the fit results from the 
quiescent data as priors was conducted, and the results 
are consistent with those when the system parameters 
were fixed at the best-fit values. 
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Fig. 6. — Disk flux contribution as a function of time. The 
squares, triangles, and diamonds indicate the g' , r 1 and i'-band 
data, respectively. 

However, we note that fine structures in the outburst 
light curves can not be well described by our model. The 
best-fit model light curves are shown in Figure [5J Com- 
paring them to the data, the second light curve, while 
sinusoidal-like, seems narrower than our model one, and 
the data points in the third and fourth ones are not ac- 
curately represented and the third one shows a second 
peak after the main peak. We conclude that modulation 
changes were detected in the outburst, but given the rel- 
atively large uncertainties in the light curves, we can not 
discriminate them from our model of an irradiated com- 
panion star by fitting. 

4. DISCUSSION 

4.1. Orbital Period Determination 

It has been shown that the orbital periodicity of 
J1808.4 can be determined from time-resolved photom- 
etry with a phas e-coherent timing technique applied 
dWang et al.l I2009T) . Usin g three observations made in 
5 days. I Wang et all ()2009[ ) have found an uncertainty of 
2.8 s on the orbital period P. With this uncertainty, 
it can be shown that such timing observations within a 
maximum time span of 3.6 month are needed in order 
not to lose the track of the optical periodicity phase. 
Because our Gemini observations in 2008 August sat- 
isfies this requirement, which further improves the ac- 
curacy of the period measurement, we tested to phase 
connect all our light curves obtained during the qui- 
escent state. The orbital period of J1808.4 has been 
found to inc rease at a rate of 1.2x l0~ 4 s yr _1 from pul- 
sar timing ( Har tman et~aT1 120081 12009ft . and this large 
value, which is not understood because it is an or- 



der of magnitude larger than that predicted by stan- 
dard theoretical calculations for such a binary, might 
be confirmed from optical observations. Given that the 
modulation was variable between different bands and 
times of observations (Figure |3]), we fit each set/band 
of the light curves individually with a sinusoidal func- 
tion m = m c + rrih sin[27r(f/P + (f>o)} (where t is time, <po 
the starting phase, and P is the orbital period fixed at 
7249.16 s), subtracted the obtained constant magnitude 
m c from them, and normalized them with the obtained 
semi-amplitudes to^. The normalized light curves were 
then fit with a single sinusoid again, and we found that 
the best-fit has orbital period P = 7249.147 ± 0.003 s, 
while x 2 = 4880 for 327 degrees of freedom. The large 
\ 2 value reflects systematic uncertainties resulting from 
photometry, intri nsic scattering of the data points from 
a sine le sinusoid C Wang et al.l l2009), errors from normal- 
izing short light curves. Adding a systematic error to the 
light curve uncertainties to lower the reduced \ 2 ~ 1, we 
found P = 7249.144 ± 0.008 and <j) = 0.675 ± 0.003 at 
MJD 54599.0 (TDB; Phase = 0.0 corresponds to the 
ascending node of the pulsar orbit), whic h are consistent 
with that obtained from pulsar timing (jHartman et al.l 
2009) within the uncertainties. Since the measurement 
accuracy of a period is generally proportional to the 
length of the time span, it is unlikely to be able to mea- 
sure the orbit change rate via long-term optical photom- 
etry of the source, unless the uncertainties of light curve 
measurements can be significantly reduced. 



model. The structure of an accretion disk in quiescence is 
known to be different from that of the standard, steady- 
state disk model (e.g..|Dubus. Hameurv. fc Lasotall2001l : 
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Fig. 7. — Broad-band spectra of J1808.4. Squares, triangles, and 
diamonds indicate the g' , r' , and %' data, respectively. The data 
points before the onset of the brightening can be fit with a disk 
model that has temperature profile T(r) oc with f = —0.5 (solid 
curve), while during the onset, £ changed to be steeper (£ ~ —0.8; 
dashed curve). A £ = —0.5 model (dotted curve) can not provide 
an acceptable fit to the data points at the onset. 

4.2. Accretion Disk During the Period 

The accretion disk components resulting from our 
MCMC fitting are shown in Figure [7] as a function of 
wavelengths. While the data points are only over 2 or 
3 bands, they definitely indicate the change of the disk 
during the period from quiescence to the onset of the 
brightening on 2008 Aug. 7. The data points in qui- 
escence can generally be fit with a simple accretion disk 



Bapti sta fc Bortolettolte004D . The temperature profile of 
the former is of the same function, T(r) oc where r 
is disk radius, but is usually flatter than that of the lat- 
ter (£ = —0.75 is the standard value for the latter case). 
Fitting the quiescent data points with the temperature 
profile, we found that T = 4300 K (at r = 10 10 cm) and 
£ = —0.5 provide the best-fit. In this fitting, parameters 
such as source distance, reddening, and inclination angle 
were fixed at the best-fit values found from our MCMC 
fitting, and the inner and outer edges of the disk were as- 
sumed to be at the light cylinder of the pulsar (~ 1.2 x 10 7 
cm) and the circularization radius (~ 3 x 10 10 cm), re- 
spectively. The spectrum of the best-fit disk model is 
shown as the solid curve in Figure [7] 

The data points on Aug. 7 are 2-4 times brighter than 
those obtained three days earlier, indicating the onset of 
the disk brightening. In addition, a spectral change from 
rising to falling is clearly visible. Although we have only 
two data point, if £ = —0.5 is fixed, the disk model can 
not provide an acceptable fit (dotted curve in Figure [3 
T — 5600 Ka.tr — 10 10 cm). Instead, we found a steeper 
profile such as £ = -0.8 (T = 5000 K at r = 10 10 cm) 
is preferred. The £ value may suggest that the disk was 
changing to the standard, steady-state case, although the 
brightness at the time (average i' ~ 19.9) was approx- 
imately 7 times low er than that at the outburst peak 
(average i' ~ 17.8; lElebert et all 120091 ). From Aug 7 
to September 22 when the outburst was first detected 
at X-rays, it is of great interest to study how the disk 
evolved during this one and a half month period. Since 
the next outbur st of J1808.4 w ill probably occur after 
2011 November (|Gallowavl l2008i). a close monitoring of 
the source near this time might help us learn the disk 
evolution over such a period. 

4.3. Modulation in the Outburst 

While our short coverage of the outburst modulation 
prevents us from conclusive analysis, we have likely de- 
tected modulation changes during the outburst. The sec- 
ond light curve is narrower than that of the normally seen 
sinusoid-like shape, probably suggesting that a smaller 
surface area of the companion star was strongly heated 
to a high temperature. Using our model, we tested this 
and found that a narrower light curve is produced but 
with a much higher 25,000 K temperature required. It is 
therefore plausible that during the outburst, a hot spot 
could have formed on the surface of the companion due 
to strong irradiation by X-rays from the central pulsar. 

The third and fourth light curves, 3 days after the sec- 
ond one, ch anged to have a flat m odulation with 0.1 mag 
amplitude. lElebert et al.l (|2009l ) have reported their de- 
tection of a slow-rise and a fast-decline modulation on 
2008 Oct. 1, which would strongly suggest the appear- 
an ce of superhump modulation since or at least on Oct. 
1 (jRetter et alj|l997t IWang fc Chakrabartvll2l)Tol ). As 
observed in cataclysmic variables (CVs), superhumps do 
often appear a few days after the start of an outburst 
(or superoutburst as in CVs). Before an accretion disk 
dumps nearly all its stored mass to the companion, the 
disk in the beginning of the outburst expands beyond 
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the 3:1 resonance radius, which induces the tidal insta- 
bility. As a result, the disk develops into an ec centric 
form , producing superhump modulation (e.g., see lOsakl 
119961 for details). However in our data, not only are the 
light curves significantly different from the Oct. 1 one 
obtained by Elebert et al. (2009; for example, its mod- 
ulation amplitude was 0.2 mag), but also no clear phase 
shift of the modulation is seen. In the last two panels of 
Figurc[2] we marked superior conjunction of the compan- 
ion (i.e., when the pulsar is at 270° from the ascending 
node) with two vertical dash-dotted lines, and as can be 
seen, the maximums of the modulation were in phase 
with superior conjunction of the companion. For super- 
hump modulation in J1808.4 and given its mass ratio 
q = 0.07/1.6, the expected phase shift would be 0.1 per 
day, estimated from the empirical relation between su- 
perhump and o rbital periods and q among known super - 
hump binaries (jPatterson et al.ll2005t IWang et alJ l2009). 

Instead, our last two outburst li ght curves are similar 
to that obtained on Oct. 7/8 by Elebert et all (|2009D . 
The origin of this 0.1 mag amplitude modulation is not 
clear. We speculate that it could be a combination of 
multiple modulations from the irradiated companion and 
the disk. Better observational coverage of the source's 
outbursts and detailed analysis of each component's pos- 
sible modulation are needed, in order to understand the 
observed low-amplitude modulation. In any case, given 
that the modulation showed significant variations over 
a few days, further observational studies of the source's 
modulation during its next outburst are warranted. 

4.4. Albedo of the Companion Star 

Day-side optical emission from the companion is due 
to reprocessing of energy flux from the central pulsar. 
The luminosities were calculated by integrating over the 
surface of the star in our model. In quiescence, they 
were in a range of (5-8) xlO 31 erg s _1 . The values are 
approximately 30% smaller than those calculated by sim- 
ply assuming a blackbody sphere for the day-side surface 
of the companion. Using the parameters of the binary 
system derived from our MCMC fitting, we can estimate 
that a fraction / = {R comp /D b ) 2 /4 ~ 0.65% of the total 
energy flux from the pulsar is received by the compan- 
ion for isotropic emission, where i? com p is the radius of 
the companion and Dt, = 6.6 x 10 10 cm is the separa- 
tion distance between the pulsar and companion. The 
known energy output from the pulsar in quiescence is its 
rotational energy loss rate L s d and its X-ray emission, 
where the latter is two orders of magnitude lower than 
the former. Therefore the energy output absorbed by the 
companion is approximately (1 — 77*)/£ s d — 6 x 10 31 erg 
s _1 , where 77, is the albedo of the companion in quies- 
cence. Comparing this value to the luminosity range of 
the companion's day-side emission and considering un- 
certainties on the binary system and the assumption of 
isotropic emission from the pulsar, 77, ~ is required. 

The optical luminosity of the companion during the 
outburst on 2008 Sept. 29 was found to be 1.4xl0 33 
erg s _1 , while the data on the other three nights are not 
considered because of the deviation of their modulations 
from a simple sinusoid-like shape. Nearly at the same 
time (MJD 54738.28), the 2.5-25 keV X-ray luminosit y 
of the pulsar was 2.5x 10 36 erg s -1 (|Hartman et aljf2 009). 



where 3.6 kpc source distance (Table [T]) is used. The X- 
ray energy received by the companion would be at least 
1.6xl0 34 erg s _1 for isotropic emission. This value is 10 
times larger than the optical luminosity, implying that 
the X-ray albedo of the companion is r]x — 0.9. Such 
a large albedo value has been found in studies of accre- 
tion disks' X-ray reprocessing (|de Jong et al.lfl996h . The 
large difference in albedo during quiescence and the out- 
burst probably provides additional evidence for having a 
different heating source in the quiescent state. 

5. SUMMARY 

We have obtained nearly simultaneous g' and i' light 
curves of J1808.4 in 2008 Aug. and r' light curves in 
2008 Sept. and Oct. during the source's outburst. In 
the former the onset of the disk brightening before the 
outburst was caught by coincidence, and in the latter 
modulation variations were seen with a sinusoidal-like 
shape changing to more-or-less a triangle-like one. 

We employed an MCMC technique to fit all the qui- 
escent data, which incl ude those preyiousl y published in 
iDelove et all (|2008| ) and lWang et"aH (|2009f ) , with an irra- 
diated companion model. We found a tight constraint on 
the inclination angle of the binary system: i = 43° ± 3° . 
The resulting constraints on the masses of the neutron 
star and companion were improved to be 1.6 35 -^0 
and 0.07to;o} M , respectively. 

From our fitting, the derived day-side temperature of 
the companion appeared to be decreasing over the pe- 
riod of ~500 days during which our quiescent data were 
obtained. We suggest that this decrease could be caused 
by a thickening accretion disk, which thus blocked an in- 
creasing part of energy output from the central pulsar. 
The derived disk components can be described by a sim- 
ple disk temperature profile, while a change to a steeper 
spectrum was seen at the onset of the disk brightening, 
indicating that the disk had already started structural 
changes 1.5 months before the outburst. 

Because of relatively poor quality of the data obtained 
during the outburst, conclusive analysis of the data can 
not be conducted. We have likely detected a narrower 
modulation peak than that seen in quiescence, which 
could be caused by the existence of a hot spot on the 
surface of the companion due to intense irradiation by 
the X-ray pulsar. The origin of the flat, low-amplitude 
modulation is not clear. In order to study it, good ob- 
servational coverage of the source during its outbursts is 
needed. 

We estimate the albedo of the companion during its 
quiescence and the outburst, and find that its values were 
approximately and 0.9 in the former and latter states, 
respectively. This large difference helps support that the 
neutron star in J1808.4 switches to be a radio pulsar in 
the quiescent state. 
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